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Angular dependent core-hole screening effects have been found in the cobalt if-edge x-ray ab- 
sorption spectrum of LiCo02, using high- resolution data and parameter- free GGA+U calculations. 
The Co Is core-hole on the absorber causes strong local attraction. The core-hole screening on the 
nearest neighbours cobalt induces a 2 eV shift in the density of states with respect to the on-site 
Is — 3d transitions, as detected in the Co K pre-edge spectrum. Our DFT+U calculations reveal 
that the off-site screening is different in the out-of-plane direction, where a 3 eV shift is visible in 
both calculations and experiment. The detailed analysis of the inclusion of the core-hole potential 
and the Hubbard parameter U shows that the core-hole is essential for the off-site screening, while 
U improves the description of the angular dependent screening effects. In the case of oxygen K- 
edge, both the core-hole potential and the Hubbard parameter improve the relative positions of the 
spectral features. 

PACS numbers: 78.70.Dm,71.15.Mb,71.20.-b 



I. INTRODUCTION 

Layered LiCo02 has been extensively studied in re- 
lation to its wide use as a cathode material in recharge- 
able Li-batteries for portable devices. Upon electrochem- 
ical cycling (Li intercalation and deintercalation) , defec- 
tive cobaltites LLjCo02 show fatigue and degradation, 
which limit the performance of the batteries. This has in- 
duced many experimental and theoretical efforts to inves- 
tigate and interpret the modifications occuring in these 
materialSfii2^i^i£i£' 7 -&2 In particular, element - and or- 
bital - selective spectroscopies, such as X-ray Absorption 
Spectroscopy, X-ray Emission Spectroscopy and X-ray 
Photoemision Spectroscopy, have been used to track lo- 
cal changes in the electronic and crystal structure, as a 
complement to X-Ray Diffraction. However, hardly any 
X-ray Absorption Near Edge Structure (XANES) data 
is available at the Co TGedge in Li :c Co02 compounds, 
and all the existing data is measured in Total Fluores- 
cence Yield, showing broad features, which are difficult to 
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Over the past years, X-ray Absorption 



Ka\ detection line yields the true absorption^ 4 Thus, 
the considerably enhanced resolution gives new opportu- 
nities to infer unique details of the electronic structure of 
Co in LiCo02, which is the goal of this paper. 

Admittedly, the fine interpretation of the absorp- 
tion features (i.e., beyond a 'finger print' analysis) re- 
quires parameter-free simulations of the spectra. First- 
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Spectroscopy using High Energy Resolution Fluorescence 
Detection (HERFD-XAS) has become a valuable tool to 
reveal subtle features of the absorption spectra, as can be 
direclty seen from Fig. Q] at the Co K-edge in LiCo02- 
HERFD-XAS is based on a two photons process (photon- 
in, photon-out) and thus may not yield true absorption 
spectra, due to strong electron-electron interactions in 
the final stated However, in the case of Co if -edge in 
LiCo02, it has been shown that HERFD-XAS using the 
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FIG. 1: Comparison between XANES spectra recorded at 
the Co Tf-edge using Total Fluorescence Yield (TFY, dashed 
line) and using High-Energy Resolution Fluorescence Detec- 
tion (HERFD-XAS, solid line). 
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principles calculations based on Density Functional The- 
ory (DFT) have already proven to be very efficient for 
many systems, both to calculate the ground-state and the 
XANES spectra. However, a first limiting factor is the 
description of the electronic interaction, which is taken 
into account in a mean way (Local Density Approxima- 
tion, labeled LDA, or General Gradient Approximation, 
labeled GGA). As an illustration, LiCo02 is known to 
be an insulator with a band gap of 2.7 ± 0.3 eV^ but 
LDA or GGA approaches underestimate the value of the 
gap in the ground state^^ To improve the description 
of the 3d-electron correlation, a solution is to include the 
Hubbard parameter Ui 2 ^ In the XANES calculation, an 
additional difficulty is to describe the core- hole effect. 
Both effects can possibly lead to a misassignement of the 
XANES features, due to wrong calculated peak positions 
and intensities. 

In this paper, we investigate the electronic struc- 
ture and the XANES of LiCo02 using a parameter- 
free DFT+U approach, which has recently been devel- 
opped. U is determined self-consistently as an intrinsic 
response of the material, and thus is no longer a fitting 
paramete n 16 i 17 This method has been already applied to 
calculate the ground state of several systems (molecules 
and crystals), showing a significant improvement to de- 
scribe the electronic structure! 17 ' 18 i 19 ' 21 i 22 XANES calcu- 
lations using DFT+U, though very promising, have only 
been performed on a restricted number of systems i 21 i 22 
Here, we use these novel theoretical approaches to in- 
terpret finely the subtle XANES features revealed by 
HERFD-XAS at the Co K-edge in LiCo0 2 . 

The paper is organized as follow. Section II gives the 
details of the XANES and of the Density Of States (DOS) 
calculations performed in LiCo02 using the GGA+U ap- 
proach. Section III is devoted to the results and the dis- 
cussion. First, we compare the ground state calculations 
obtained by GGA and GGA+U (Sec. IIIA). Then, after 
presenting the XANES spectra calculated at the Co K- 
edge, we interpret the spectral features using local DOS 
and interpret the effect of U and of the core-hole on the 
theoretical spectrum (Sec. IIIB). In Sec. IIIC, we focus 
on the Co if pre-edge region and interpret the origin of 
the local and non-local features. Finally, for complete- 
ness, we compare the theoretical XANES spectrum at 
the O if -edge to experimental data recorded in Total 
Fluorescence Yield (Sec. HID). Conclusions are given in 
Section IV. 



II. COMPUTATIONAL DETAILS 

First principles calculations were performed using the 
Quantum-Espresso first-principles total-energy code* 2 ^ 
The code uses plane waves and periodic boundary con- 
ditions. The XANES spectra are obtained in two steps: 
first the charge-density is obtained self-consistently using 
the PW package of the Quantum- Espresso distribution, 23 
then the XANES spectrum is computed in a continued 



fraction approach using the XSPECTRA packagej 22 ' 24 i 25 
As Co and O require large cutoffs to be simulated 
with standard norm-conserving pseudopotentials, we 
use ultrasoft pseudopotentials with two projectors per 
channel* 2 -^ The electronic configurations are 3d 8 4s 1 with 
non-linear core correction 27 for Co, 2s 2 2p 4 for O and 
2s 1 2p° with non-linear core correction for Li. The use 
of Ultrasoft pseudopotentials limited the cutoff of the 
plane-wave expansion at 45 Ry energy. The generalized 
gradient approximation (GGA) was adopted. 28 

When no core-hole was considered in the final state, 
a 12-atoms hexagonal unit cell with experimental lattice 
parameter and atomic positions was usedi 2 ^ Electronic 
integration in the self-consistent run was performed on a 
8x8x8 Monkhorst-Pack fc-point grid shifted along the c- 
axis. The DOS and Lowdin projected-DOS calculations 
were performed using al2xl2xl2 fc-points mesh and 
a Gaussian broadening of 0.3 eV. 

Core-hole effects were treated in a supercell approach 
using a 3x3x1 hexagonal supercell containing 108 atoms 
(54 O, 27 Li, 26 Co and one absorbing Co) including the 
core-hole in the absorbing- atom pseudopotential. The 
distance between two neighboring absorbing atoms, upon 
the application of periodic boundary conditions, is 8.5 A, 
which is large enough to avoid interactions. In this case 
electronic integration in the self-consistent run was per- 
formed using a 4 x 4 x 4 fc-point Monkhorst-Pack mesh 
shifted along the c-axis. The XANES calculation with a 
core-hole in the final states required electronic integra- 
tion over a 4 x 4 x 4 fc-point grid. 

In the study of the electronic properties of LiCo02 the 
use of a Hubbard U parameter was found essential in 
order to reproduce the XAS spectrum and the magni- 
tude of the electronic gap. In particular we impose U 
on the d-states of the Co atoms and calculate its value 
self-consistently (U sc f) using the method of Refs 
We obtain U sc f=5.6 eV. Thus there are no free parame- 
ters in our calculation. The XANES Co and O if -edges 
cross-sections, both in their electric dipole (El ) and elec- 
tric quadrupole parts (E2) were then calculated includ- 
ing the Hubbard U termini The isotropic spectra were 
calculated in the correct point group symmetry (3m for 
Co, 3m for O) according to the procedure given in Ap- 
pendix. At the Co if -edge, a constant broadening pa- 
rameter (0.55 eV) was used in the continued fraction. 
This value is smaller than the broadening due to the Co 
Is core-hole lifetime: indeed, the calculation is compared 
to experimental HERFD-XAS data recorded by monitor- 
ing the absorption using the Ka\ fluorescence line, which 
allows to remove partly the core-hole lifetime!^ 8 . The ex- 
perimental data is taken from Ref. [3 For the O if -edge, 
as the experimental data was measured in Total Fluores- 
cence Yield^ we used a broadening of 0.6 eV between 
525 eV and 532 eV, and a broadening of 0.8 eV above 
532 eV. 
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FIG. 2: Density of states of LiCoCh calculated in GGA (pan- 
nel (a)) and GGA+U (pannel (b)), together with Lowdin pro- 
jected density of states on Li, Co and O atoms. The Fermi 
level is set at the bottom of the conduction band. 



III. RESULTS AND DISCUSSION 

A. Ground state calculations 

The total and Lowdin projected DOS have been cal- 
culated in GGA and GGA+U for the ground-state of 
LiCo02 (i.e., with no core- hole). They are shown in 
Fig. [2] (pannels (a) and (b), respectively). We used 
spin-polarization and started from an initial configura- 
tion, where each Co atom carries either a magnetic mo- 
ment equal to zero, or slightly positive (0.2 /ib)- In both 
cases, we observed a convergence to a non-magnetic state, 
where the final magnetic moment on the Co atoms is 
zero. This is consistent with previous experiments and 
calculations, which show that Co 3+ in LiCo02 is in a low 



spin state S = 0&££L For U = 0, DFT converges to an 
insulating ground state with a 1.1 eV electronic gap to 
be compared with the experimental one 2.7 ± 0.3 eV£ 
The Lowdin projected DOS show that the valence band 
is dominated by Co and O states, while the conduction 
band is dominated by Li and Co states. The gap is be- 
tween t2 g and e g Co states slightly hybridized with O 2p 
states. These results are consistent with previous DFT 
calculations on LiCo02^ i ^ 

DFT+U calculations using our calculated Hubbard 
parameter lead to a 2.3 eV gap (see panel (b) in 
Fig. [2]), in good agreement with the one determined 
from x-ray photoemission and bremsstrahlung isochro- 
mat spectroscopy. 5 U shifts to lower energies the occu- 
pied states close to the Fermi level, i.e., Co 3d states plus 
Li p and O p states. For the empty states, the position 
in energy with respect to the Fermi level is identical to 
the calculation with U = 0. 



B. XANES at the Co K-edge 

The calculated and experimental isotropic Co i*T-edge 
spectra for LiCo02 are shown in Fig. [3l The experimen- 
tal HERFD-XAS spectrum (thin solid line) shows a set of 
three features at 7717.8 eV (peak a), 7720.9 eV (peak b) 
and 7722.2 eV (peak c), and the main edge is composed 
of a broad feature (peak d) centred at 7726.9 eV. This 
spectrum is compared to theoretical spectra obtained us- 
ing four different GGA calculations, i.e., by switching 
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FIG. 3: Comparison in the edge region between experimental 
(thin solid line) and theoretical isotropic XANES spectra at 
the Co -ftT-edge in LiCo02, obtained by four different calcu- 
lations: GGA+U with or without the Co Is core-hole (black 
and orange solid lines), GGA with or without the Co Is core- 
hole (black and orange dashed lines). 
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on/off the Hubbard parameter U and on/off the Co Is 
core-hole. These four spectra have been shifted in energy, 
so that the position of the main edge (maximum of the 
absorption) matches with the experimental one. 

First, all the experimental features are reproduced in 
the calculations. As mentioned in the introduction, com- 
puting the absorption cross-section is here fully justified, 
since at the Co if -edge in LiCo0 2 , the HERFD-XAS 
spectrum is a measurement of the true absorption^ By 
monitoring the influence of U and of the core-hole in 
the calculation, we can therefore draw conclusions on the 
core-hole effects in the material investigated experimen- 
tally. 

Second, when comparing the four theoretical spectra, 
differences in the peak positions and relative intensities 
are present. The best calculation is obtained when tak- 
ing into account both the Hubbard parameter on the Co 
3d orbitals and the Co Is core- hole (black solid line). 
Indeed, the relative intensities and positions of peaks b 
and c, as well as the shape of peak d, are particularly im- 
proved. The Hubbard parameter has no effect when the 
core-hole is not included (orange and black dashed lines), 
since the relative positions of the peaks above the Fermi 
level are identical (see Fig. [2]). The main effect of the 
core-hole for U = is to improve the shape of peak d, by 
merging the two subfeatures into a sharper one. When 
the core-hole is on, the Hubbard parameter mainly af- 
fects the intensity of peak a, the relative intensities and 
positions of peaks b and c, leading to a better agreement 
with experiment. 

In order to interpret the influence of U, of the core-hole 
and their possible anisotropic effects, we have calculated 
the electric dipole cross-sections for a polarization vector 
either perpendicular or parallel to the c axis of the crystal 
(Figs. 2k andQ})). Indeed, as the spectra shown in Fig. [3] 
are isotropic, they show an average absorption and thus 
they may be misleading. The dichroic XANES spectra (e 
perpendicular and e parallel) correspond respectively to 
transitions to the empty (p x , p y )- andp z -states projected 
on the absorbing Co atom. Although no experimental 
data is available, the theoretical spectra enable a better 
understanding. In order to compare them, the spectra 
have been shifted in energy to match the position of the 
feature at 7730 eV (Fig. S^). This feature corresponds 
indeed to delocalized Co p states, which we assume to 
be less affected by U and the core-hole. U has no effect 
on the isotropic theoretical spectrum when the core-hole 
is not included. This can be explained by the fact that 
U is added on Co 3d orbitals, which are very localized 
states and do not contribute to the edge. When the Co Is 
core-hole is added, the intensity of peak d is increased for 
both polarizations, together with a shift to lower energies. 
Both effects can be explained by the core-hole attraction 
and it affects in a similar way p x , p y and p z states. This 
accounts for the fact that peak d becomes narrower in the 
isotropic spectrum, leading to a better agreement with 
experiment. When both the Hubbard parameter and the 
core-hole are included, peaks a, b and c are significantly 
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FIG. 4: Electric dipole cross-section at the Co K-edge calcu- 
lated using the four different approaches. Pannels (a) and (b) 
show the cross-sections calculated for an incident polarization 
vector perpendicular and parallel to the c axis. 



affected. Peak a is shifted up and peak c is split into 
two subfeatures, which consequently makes peak b grow. 
As this was not present for U = with the core-hole, 
it can be attributed to the Hubbard parameter. These 
effects are similar for both polarizations. However, this is 
different for peak d, which is shifted to low energies and 
which is more intense, but only for the corresponding 
p z states. From the previous analysis, we can conclude 
that the Co Is core-hole mainly affects the intensity of 
the XANES features by attracting Co empty p states to 
lower energy. The Hubbard parameter on Co 3c? states 
affects the peak positions, but only when the core-hole is 
present. This can be explained by the fact that the core- 
hole attraction mixes Co and O empty states in such a 
way that U can act also on the Co p empty states. 

Finally, we want to comment on the fact that our DFT 
calculations reproduce all the experimental XANES fea- 
tures. Indeed, DFT can only take into account one sin- 
gle Slater determinant and thus neglects charge-transfer 
effects. When the latter are significant, additional fea- 
tures are visible in the edge of the experimental spec- 
trum (e.g., Cu K-edge in La2Cu04)22, but they cannot 
be calculated using a single-particle approach. The ex- 
istence of charge transfer induced transitions can be in- 
ferred from Is X-ray Photoemission Spectroscopy (XPS) 
spectra. If intense satellite features are visible, one ex- 
pects a modification of the XANES K-edge features. The 
good agreement obtained between theoretical and exper- 
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FIG. 5: Theoretical Co Is XPS cross-section calculated using 
the charge-transfer multiplet approach. 



imental XANES spectra in LiCoC>2 enables already to 
conclude that charge transfer satellite transitions are neg- 
ligible. To support this result, we have calculated the Co 
Is XPS cross-section using the charge-transfer multiplet 
theory. This calculation used the values of the crystal 
field, the Slater integrals and the charge transfer param- 
eters from Ref. |5| It included two configurations in the 
initial and final states of XPS: 3d 6 and 3d 7 L, where L 
denotes a hole on an oxygen ligand. The theoretical Co 
Is XPS cross-section is shown in Fig. [5] We found that 
95 % of the intensity goes to the main peak, which is 
visible at 7709 eV in binding energy and which is the so- 
called well-screened peak, a mixture of mainly Is 3d L 
with some ls x 3(i 6 final states. As the high-binding en- 
ergy satellite features are very weak, satellite transitions 
in the XANES will not be visible, which is in line with 
our results. 



C. Insight from the Co K pre-edge 

1. Comparison between experiment and calculations 

We now do a similar comparative analysis in the Co 
K pre-edge region (shown in Fig. [6]). The experimen- 
tal isotropic pre-edge spectra (thin solid line) shows two 
main features, labeled a and /3, which are visible respec- 
tively at 7708.9 eV and 7710.9 eV. Additionnally, a small 
shoulder labeled 7 is visible at 7712.3 eV. At this stage 
of the paper, we will focus the analysis on the intense 
peaks a and (3, since peak 7 is not discriminative of the 
different calculations (the broadening used in the calcu- 
lation being too large). The origin of this shoulder will 
be discussed later in the paper. 

All the theoretical spectra are shifted to high energy 
compared to experiment (the position of the edge being 
the same). Such a shift is systematically observed at the 
K pre-edge of transition metal ions i 21 i 33 ' 34 i 35 it is due 
to the fact that the Is core-hole is overscreened in the 
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FIG. 6: Comparison in the Co K pre-edge region between ex- 
perimental (thin solid line) and theoretical isotropic XANES 
spectra in LiCo02, obtained by four different calculations: 
GGA+U with or without the Co Is core-hole (black and or- 
ange solid lines), GGA with or without the Co Is core- hole 
(black and orange dashed lines). 



calculation, which consequently underestimates the split- 
ting between the electric dipole and electric quadrupole 
transitions. Nevertheless, as this shift is observed for all 
the theoretical spectra, we will now compare the relative 
intensities and peak positions. There are significant dif- 
ferences between the four theoretical spectra. When the 
Hubbard parameter and the core-hole are not included 
in the calculation (orange dashed line) , the pre-edge con- 
sists of one single broad peak at 7711.1 eV. When U is 
added (orange dashed line), this feature is shifted only 
by 0.1 eV to lower energies. If the core- hole is taken into 
account, either with U = or with U = U sc f (orange and 
black solid lines) , this single broad intense feature is split 
into two sharp ones with lower intensity. The splitting 
is higher for U = U sc f, which yields the best agreement 
obtained between calculation and experiment. We can 
already conclude that for the Co K pre-edge region, the 
core-hole affects the splitting of the two main peaks, and 
U influences strongly the splitting of the two, as well as 
their positions with respect to the edge. For U = U sc f 
and the core-hole, the intensity of peak a is however un- 
derestimated by 40 % with respect to experiment, while 
that of peak f3 is well reproduced. The splitting between 
peaks a and j3 is also smaller than in the experimental 
data (1.4 eV vs 2.0 eV). We will further comment these 
points after analyzing in detail the origin of the pre-edge 
features. 

The origin of the pre-edge features can be clarified by 
plotting separately the electric dipole (El ) and the elec- 
tric quadrupole (E2) contributions. Figure^ shows the 
El (solid line) and E2 (dashed line) contributions, as 
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it improves visibly the intensity ratio between the peaks. 
Therefore, only the relative intensities should be com- 
pared to the experimental ones. From Fig. [7Ji, we can 
conclude that peak a is due to pure electric quadrupole 
transitions (Is — 3d) and that peak (3 is due to almost 
pure electric dipole transitions (Is — p). Peak 7 is due to 
pure electric dipole transitions. 

To go further into the analysis, we have calculated the 
angular dependence of the electric dipole and the elec- 
tric quadrupole cross-sections. This can be achieved by 
calculating the absorption cross-section for various orien- 
tations of the polarization and wave vectors, with respect 
to the crystal. Figure [7b shows that peak a corresponds 
to 3d states that have, non-surprisingly, a pure e g charac- 
ter. Indeed, when the transition operator £ ■ rk • r is cho- 
sen along a Co-0 bond (i.e., e g -like orbitals are probed), 
the intensity of peak a is maximized. When it is cho- 
sen between two Co-0 bonds (i.e., ^g-like orbitals are 
probed), peak a disappears almost completely (Fig. LTJd). 
Such X-ray Linear Natural Dichroism effects are already 
well-known and are used to derive local information on 
the electronic and crystallographic structured However, 
in the case of Co 3+ , they are particularly emphasized, 
since the e g orbitals are the only ones to be empty. Be- 
cause the Co site point group symmetry is not Oh but 
slightly lower (D34), it is not possible, in principle, to ex- 
tinct completely peak a because of a slight hybridization 
between e g and ti g states. However, we found out that 
peak a has almost a pure e g character, which shows that 
the distortion of the Co site from perfect octahedral sym- 
metry is very small. Figure [7J; shows that peak has a 
mixed p x , p y , p z character and the p z states are at slightly 
lower energy. Peak 7 has a dominant (p x ,p y ) character: 
its intensity can thus be enhanced by doing measurement 
on a single crystal using a polarization vector perpendic- 
ular to the c axis. 



FIG. 7: (a) Comparison between experimental (thin solid line) 
and theoretical isotropic pre-edge spectra (solid line) at the 
Co K-edge in LiCo02- For the theoretical spectrum we show 
the electric dipole contribution (El, dashed line), the elec- 
tric quadrupole contribution (E2, solid line) and their sum 
(E1+E2, orange solid line), (b) Angular dependence of E2 
cross-section, (c) Angular dependence of El cross-section. 
The theoretical spectrum has been shifted to the experimen- 
tal data in order to match the position of peak a, for a more 
direct comparison of the peak position. 



well as the sum of the two (orange solid line), for the cal- 
culation performed in GGA+U and with the core-hole. 
The theoretical spectrum has been shifted to the experi- 
mental data in order to match the position of peak a, for 
a more direct comparison of the peak position. Note that 
in this case, the calculation has been performed with a 
smaller broadening parameter (0.2 eV) than in Fig. [SJ in 
order to separate clearly the three peaks. In particular, 



2. Assignement of the pre-edge features 

To assign these transitions, the DOS for the supercell 
including the Co Is core-hole and U = U sc f has been 
projected on the absorbing Co atom and selected neigh- 
boring atoms (O, Li and Co) in a range of 5 A from the 
absorber. The local DOSs are plotted in Fig[8l We dis- 
tinguish: the first and second oxygen neighbors (O in 01 
plane at a distance of 1.92 A and 3.41 A respectively), 
the first and second Co neighbors in the same plane as 
the absorber (Co in plane, at a distance of 2.82 A and 
4.88 A respectively), the first Li neighbors (Li plane, at 
a distance of 2.85 A), the first Co neighbors in the plane 
below the absorber (Co out of plane, at a distance of 
4.96 A), and the first O neighbors between the Li plane 
and the Co out of plane (O in 02 plane, at a distance of 
3.66 A). 

As can be seen from the local d-DOS on the absorb- 
ing Co (Fig. [5^, solid line), peak a is due to electric 
quadrupole transitions from the Is orbital to 3d states lo- 
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calized on the Co absorber (intrasite excitations). More- 
over, although empty 2p states of the O nearest neighbors 
are found at the same energy position (Fig. [8^, dashed 
line), we found no evidence of direct transitions to these 
states, contrary to the case of Ni K pre-edge in NiOi^i 

The local p-DOS on the absorbing Co atom (Fig. 
dashed line) shows that peak (3 is mainly due to electric 
dipole transitions, i.e. from the Is orbital to Apt empty 
states. At this energy (2 eV above the Fermi level), there 
are actually two different contributions slightly shifted 
one from each other: one centred at 1.6 eV (peak (3') and 
a second one, centred at 2.0 eV (peak j3"). This accounts 
for the fact that peak (3 is broader than peak a. The local 
DOS on the neighbors lying in the same plane as the ab- 
sorber clarify the origin of these two subfeatures. Peak (3' 
is due to transitions to empty Ap states of the absorber, 
hybridized with 2p states of the O nearest neighbors (Fig. 
IHfe, dashed line) and 3d states of the Co nearest neigh- 
bors (Fig. [SJs, solid line). Peak f3" is due to transitions to 



empty Ap states of the absorber, hybridized mainly with 
2p states of the O second-nearest neighbors (Fig. [5f, 
dashed line) and with 3c? states of the Co second near- 
est neighbors (Fig. [5J;, solid line). The splitting between 
peaks f3' and (3" is due to a different core-hole screen- 
ing between the (Co,0) first neighbors and the (Co,0) 
second neighbors: as the latter lie further from the ab- 
sorber, the Co Is core-hole is more screened. Therefore, 
peak (3 has a double non-local [off-site) origin, which is 
revealed by the core-hole attraction. We have found that 
peaks j3' and j3" both owe their broadening to the fact 
that the Co Ap z empty states lie at lower energies than 
the (Ap x , Ap y ) states, as shown in Fig. [7];. However, this 
is not due to a different screening of the core-hole along 
these directions, since in the ground-state (i.e., with no 
core-hole) , the Co p z empty states already lie lower in en- 
ergy (0.2 eV) than the empty p x ,P y states. Therefore, it 
is the strong anisotropy of the structure that introduces 
direction-dependent hybridizations and lifts orbital de- 
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generacy, even when the core- hole is absent. 

Peak 7 is due to non-local excitations, i.e., to on-site 
Co empty 4p states hybridized with empty 3d states of 
the Co atoms, which are in the plane below and above 
the absorber (Co out of plane, Fig. [Hi). Such far off-site 
transitions (at 2.6 eV above the Fermi level) are medi- 
ated by the Li 2p empty orbitals (Fig. [7fr, dashed line) 
and the 2p states of the O atoms (nearest neighbors plane 
02, dashed line), which are between the Li and the Co 
out of plane atoms. These Co atoms lie a distance of 
4.96 A from the absorber, which is close to the distance 
between the absorber and the Co second neighbors in 
plane. This implies that the distance being the same, the 
core-hole is more screened along the c direction than in 
the perpendicular plane. As a consequence, the splitting 
between peaks f3 and 7 is due to an angular dependent 
screening of the core-hole. This original result can cer- 
tainly be attributed due to the anisotropic layered struc- 
ture of LiCoC>2. As we mentioned before, the far off-site 
peak has mainly a (p x , p y ) character. It can appear at 
first sight surprising that out of plane hybridization is 
achieved mainly by the O in plane p-orbitals, taking into 
account the strong anisotropy of the structure. However, 
a look at the structure reveals that the c (or z) axis is 
directed along the local C3 axis of the C0O6 octahedra. 
This implies that the Co-0 bonds are not directed along 
the x, y, and z directions and that the p x , p y or p z char- 
acter discussed does not correspond to the traditionnal 
reference frame. In such a description, it is thus not sur- 
prising to have a p Xl p y character in the peak probing the 
Co atoms out of plane. 



3. Influence of U and of the core-hole 

As mentioned before, both the Hubbard parameter and 
the core-hole must be included in the calculation of the 
Co if pre-edge spectrum, in order to obtain a good agree- 
ment with experiment. The different features are well 
separated and this is due to a complex combined effect. 
When the core-hole is off, we could expect that U shifts 
the empty Co 3d states up towards the edge, since the 
latter consists of Co p states not affected by U. This is 
however not the case. When the core-hole is switched on 
for U = 0, peaks a, (3 and 7 are separated by the core-hole 
attraction. This different behavior is due to the non-local 
character of peaks (3 and 7. The splitting between a and 
(3 increases, since the screening is all the more efficient 
that the Co and O neighbors are far from the absorber. 
The energy difference between [3 and 7 increases as well, 
because core-hole effects show an angular dependence (in 
plane and out of plane). 

When the core- hole is present, the Hubbard parameter 
on Co 3d orbitals increases again the energy difference 
between the in-site and the off-site peaks, on one hand, 
and between the two off-site peaks, on the other hand. 
This is similar to the case of Ni if -edge in NiOi^I U 
has an effect since the 3d empty states of the absorber 



and of its neighbors have already been separated by the 
differential core-hole screening. This is in line with the 
fact that U has hardly no effect if the core-hole is absent. 

The fact that the intensity of peak a is slightly un- 
derestimated compared to the two other ones raises the 
question of a possible missing electric dipole contribution 
at this energy. This can have two possible origins: (i) the 
quality of the sample investigated experimentally, (ii) vi- 
brational effects - missing in the calculation-, which intro- 
duce an electric dipole contribution in peak a. First, the 
sample is finely grinded. As most powders remain tex- 
tured and the crystal structure is trigonal, there could be 
angular dependence effects. However, these effects should 
be rather small, and we do not think this could explain 
the observed disagreement. We recall that the calculated 
spectrum is the isotropic one, and thus should be directly 
compared to the experiment. Second, another possibility 
is the vibrational effects, since the experimental data was 
measured at room temperature. At T = OK, the fact 
that the Co site is centrosymmetric (D 3 d point group) 
prevents any static mixing between the Co 3d states and 
the O 2p states in the calculated XANES spectrum. How- 
ever, lattice vibrations can remove this rule and thus, 
significant changes in the XANES spectrum may be pos- 
sible. In the case of Co if -edge in LiCo02, vibrations 
could introduce an electric dipole contribution in peak 
a. To verify the validity of this conjecture, HERFD- 
XAS measurements at the Co if -edge as a function of 
temperature should be performed and combined with a 
consequent analysis including the vibrational effects^ 



D. XANES at the O If -edge 

For completeness, we have calculated in DFT+U the 
XANES at the O if -edge, although the experimental 
data has been measured in Total Fluorescence Yield. 5 
The experimental isotropic O if -edge spectrum is shown 
in Fig. [9] (solid line with filled circles). It shows a well- 
defined feature at 528.3 eV (peak A) and two broad ones 
at 536.4 eV (peak B) and 539.9 eV (peak C). It is com- 
pared with the theoretical ones calculated in the electric 
dipole approximation (i) in GGA+U with or without the 
O Is core-hole (black and orange solid lines, respectively) 
and (ii) in GGA with or without the O Is core-hole (black 
and orange dashed lines, respectively). The four theoret- 
ical spectra have been shifted in energy in order to match 
the position of peak A with the experimental data. 

The best agreement has been obtained for the calcula- 
tion using U = U sc f and the core-hole (black solid line): 
although the intensity of peak A is underestimated, the 
relative position of the peaks is improved compare to the 
calculations without the core-hole or without the Hub- 
bard parameter. This is also the case for the shape of 
peak C and the intensity ratio between peaks B and C. 

The main effect of the core-hole is to decrease by a fac- 
tor of two the intensity of peak A, which corresponds to 
transitions to empty O 2p states. This is due to the fact 
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that part of these states are shifted down to -10 eV below 
the Fermi level (not shown), due to a strong underesti- 
mation of the core-hole screening. Similar overestimation 
of core-hole effects has also been reported in the calcula- 
tions performed at the K-edge of light elements such as 
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B and ) 37 i 39 where it was found that a half core-hole 
provides a better agreement in the region just above the 
edge. Such a drastic influence has also been shown at the 
Cu K-edge^, but only for the electric quadupole cross- 
section. Here, we show that the core-hole effects are not 
well described also for the electric dipole cross-section, 
since the intensity of peak A is better reproduced with- 
out the core-hole (though slightly overestimated). The 
core-hole improves the shape of peak C and the intensity 
ratio between peaks B and C, which can be understood 
by plotting the electric dipole cross-sections for a polar- 
ization vector e either perpendicular or parallel to the c 
axis (Figs. [T0k and [TUb) . The attraction of the core- 
hole induces a shift by -1.2 eV of the (p x , p y ) states at 
540 eV, and an overall decrease in intensity of the corre- 
sponding p states. Part of the p z states from peak C have 
been transfered to peak B (536 eV). For p x , p y states, it 
is likely that the core-hole shifts them below the Fermi 
level, since peak B has the same intensity with and with- 
out the core-hole and the intensity of peak C decreases. 

The Hubbard parameter U has no effect on the theo- 
retical spectrum when the core-hole is absent, which is 
in line with the empty DOS shown in Figf^l However, it 
does have an effect when the core-hole is on and affects 
the peak positions, similarly to the case of Co K-edge. 
Although U is added on the Co 3d orbitals, it has an 
effect on the empty O p states because they are slightly 
hybridized (Fig|2]) and because the core-hole introduces 
differential attraction on the O empty states. This results 
in a shift to higher energies of peaks B and C relatively 
to peak A. 



IV. CONCLUSIONS 

For the first time, the XANES spectra at the Co K- 
edge and O if -edge have been calculated in LiCo02 using 
a parameter-free GGA+U approach. A good agreement 
has been obtained with the experimental data, which has 
enabled to draw conlusions of several types. 

First, from the computational point of view, we have 
shown that the Hubbard parameter U improves the the- 
oretical spectra in the edge region, for both the relative 
intensities and positions of the spectral features. This im- 
provement is even more striking at the Co K pre-edge, 
where both the core-hole and U are needed to interpret 
the spectral features. The core-hole effects were found to 
be the factor preventing a better agreement with exper- 
iment: in the calculation, where we considered a static 
core-hole, the O Is core-hole was found to be too at- 
tractive, while the screening of the Co Is core-hole is 
overestimated. A better treatment of the core-hole ef- 
fects could be achieved using the Bethe-Salpeter equa- 
tion, which treats electron and hole dynamics ab initio, 
as well as electron- hole interactions^ 

Second, we have shown that the GGA+U approach 
improves significantly the description of the electronic 
structure for the ground state in LiCo02, since the value 
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of the theoretical electronic band gap is, for the first time, 
consistent with the experimental one. This implies that 
this approach could yield better results on the electronic 
structure of Li-defective cobaltites LLjCoC^. 

Third, high-resolution XANES data at the Co K pre- 
edge has revealed new subtle features, which have been 
interpreted thanks to DFT+U calculations. In the pre- 
edge, in addition to classical on-site ls-3d transitions, we 
have shown the existence of two different Co Ap - Co 3d 
intersite hybridizations, in plane and out of plane. This 
is achieved via a strong Co 3d - O 2p hybridization in the 
Co planes, and via a strong Co 3d - O 2p - Li 2p hybridiza- 
tion in the direction perpendicular to the Co planes. Al- 
though it can be argued that the presence of the core- 
hole prevents to draw any conclusion on the ground-state 
properties, we want to point out that these pre-edge fea- 
tures can only be resolved because the Co Is core- hole 
is present: the screening of the core-hole is indeed de- 
pendent both on the distance and on the direction. This 
means that HERFD-XAS measurements performed on 
a single crystal, coupled to GGA+U parameter-free cal- 



culations, can be used to track fine changes in orbital 
hybridization, such as the ones induced by electrochem- 
ical cycling in the electronic and crystallographic struc- 
ture of Lia;Co02 compounds. For example, as the in- 
crease in Li- vacancies leads to the augmentation of the c 
parameter it should affect the relative intensities and 
peak positions of the Co pre-edge features. 
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V. APPENDIX 

A. Calculation of the isotropic electric dipole 
cross-section 

For a general symmetry, the isotropic electric dipole 
cross-section is given by (a xx + a yy + a zz )/3, where <j xx , 
Gy V and a zz are the electric dipole absorption spectra 
calculated along three perpendicular directions^ 2 . Wc 
now consider the sites of oxygen and cobalt. 

In the space group Rim described in the hexagonal 
setting, there are two equivalent oxygen atoms in position 
6c: (0 z) and (0 -z). The four other O sites can be 
obtained by applying the lattice translations (1/3, 2/3, 
2/3) and (2/3, 1/3, 1/3). The point group symmetry is 
3m (or C^v). Therefore, we take z along the c axis of the 
crystal and the local symmetry group gives us a xx — a yy . 

The Co atom is in position 6c: (0 1/2). The point 
group symmetry is -3m (or D^d)- Thus, we also have 
O xx — ^yy~ 



B. Calculation of the isotropic electric quadrupole 
cross-section at the Co if -edge 

The space group symmetry operations of the group 
R3m that leave the Co site (0,0,1/2) invariant arc 

(x,y,z), (-y,x-y,z), (-x + y,-x,z), (y,x,-z), (x - 
y,-y,-z), (-x, -x + y, -z), (-x,-y,-z), (y, -x + 
V,-z), (x - y,x,-z), (-y,-x,z), {-x + y,y,z) and 
(x,x- y,z). 

The general angular dependence of electric quadrupole 
transitions in terms of a second-rank spherical tensor 
0-^(2, m) and a fourth-rank spherical tensor a^(A,m) was 
given in ref. H3- The trigonal symmetry of the Co site 
implies that the only non-zero component of er^(2, m) is 
0^(2,0), which is real. To determine the non-zero com- 
ponents of a® (A, m) for the Co site, we apply the formula 
derived in ref. |43f 

<r«(4,m) = ±J2 T ™' D i'™(Rl 

Rm' 

where R runs over the rotation part of the 12 operations 
listed above and where T^, is an Hermitian fourth-rank 
tensor (without symmetry). This gives us the only non- 
zero components of c < 5(4, m), which are cr^(4, 0) (a real 
number) and cr'^(4, ±3) (a purely imaginary number). 
We define the real a Qi {4, 3) by a Q (4, ±3) = ia Qi (A, 3). 

In a reference frame where the x axis (the z-axis, re- 
spectively) is along the a vector (the c vector, respec- 
tively) of the hexagonal lattice, this gives us the angular 
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dependence 

a Q = a Q (0,0) + -y/J(3sin 2 6»sin 2 ^- l)<x Q (2,0) 

- (35 sin 2 cos 2 cos 2 -0 + 5 sin 2 sin 2 ^ - 4) 
xcr <5 (4,0) 

— vlO sin 6 [(3 cos 2 0—1) sin ip cos cos 30 
+ cos (9(2 cos 2 (9 cos 2 r/> - 1) sin 30] cr Qi (4, 3) . 

This expression is similar but not identical with the one 
given in ref. |42| for the symmetry D^d- The difference 
comes from the fact that the relation between the Carte- 
sian axes and the local symmetry axes is different in the 
two cases. 

The isotropic quadrupolc absorption is cr** (0,0). 
In general, we need to calculate four different ori- 
entations to deduce the isotropic quadrupole spec- 



trum. However, we could obtain it from only 
three directions: ei = (—l/y/6,—l/y/2,l/\/3),kx = 
(l/VB,-l/V5,-l/>/3), £2 = (0,l,0),k 2 = (-1,0,0) 
and e 3 = (0, l/y/2, -1/V2), k 3 = (0,1/72,1/^2). The 
corresponding absorption cross-sections are 



^(0,0) 



<r«(4,0), 



a Q (0, 0) + -^a Q (4, 0) + 2\jJ-a Q (2 



14 V 14 
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The isotropic quadrupole spectrum is 
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